Linkage disequilibrium mapping in isolated populations: The example of Finland revisited

Linkage disequilibrium analysis can provide high resolution in the mapping of disease genes because it incorporates information on recombinations that have occurred during the entire period from the mutational event to the present. A circumstance particularly favorable for high-resolution mapping is when a single founding mutation segregates in an isolated population. We review here the population structure of Finland in which a small founder population some 100 generations ago has expanded into 5.1 million people today. Among the 30-odd autosomal recessive disorders that are more prevalent in Finland than elsewhere, several appear to have segregated for this entire period in the “panmictic” southern Finnish population. Linkage disequilibrium analysis has allowed precise mapping and determination of genetic distances at the 0.1-cM level in several of these disorders. Estimates of genetic distance have proven accurate, but previous calculations of the confidence intervals were too small because sampling variation was ignored. In the north and east of Finland the population can be viewed as having been “founded” only after 1500. Disease mutations that have undergone such a founding bottleneck only 20 or so generations ago exhibit linkage disequilibrium and haplotype sharing over long genetic distances (5–15 cM). These features have been successfully exploited in the mapping and cloning of many genes. We review the statistical issues of fine mapping by linkage disequilibrium and suggest that improved methodologies may be necessary to map diseases of complex etiology that may have arisen from multiple founding mutations.

Evolution of the Friedreich’s ataxia trinucleotide repeat expansion: Founder effect and premutations

Friedreich’s ataxia, the most frequent inherited ataxia, is caused, in the vast majority of cases, by large GAA repeat expansions in the first intron of the frataxin gene. The normal sequence corresponds to a moderately polymorphic trinucleotide repeat with bimodal size distribution. Small normal alleles have approximately eight to nine repeats whereas a more heterogeneous mode of large normal alleles ranges from 16 to 34 GAA. The latter class accounts for ≈17% of normal alleles. To identify the origin of the expansion mutation, we analyzed linkage disequilibrium between expansion mutations or normal alleles and a haplotype of five polymorphic markers within or close to the frataxin gene; 51% of the expansions were associated with a single haplotype, and the other expansions were associated with haplotypes that could be related to the major one by mutation at a polymorphic marker or by ancient recombination. Of interest, the major haplotype associated with expansion is also the major haplotype associated with the larger alleles in the normal size range and was almost never found associated with the smaller normal alleles. The results indicate that most if not all large normal alleles derive from a single founder chromosome and that they represent a reservoir for larger expansion events, possibly through “premutation” intermediates. Indeed, we found two such alleles (42 and 60 GAA) that underwent cataclysmic expansion to pathological range in a single generation. This stepwise evolution to large trinucleotide expansions already was suggested for myotonic dystrophy and fragile X syndrome and may relate to a common mutational mechanism, despite sequence motif differences.

A complete genome screen for genes predisposing to severe bipolar disorder in two Costa Rican pedigrees

Bipolar mood disorder (BP) is a debilitating syndrome characterized by episodes of mania and depression. We designed a multistage study to detect all major loci predisposing to severe BP (termed BP-I) in two pedigrees drawn from the Central Valley of Costa Rica, where the population is largely descended from a few founders in the 16th–18th centuries. We considered only individuals with BP-I as affected and screened the genome for linkage with 473 microsatellite markers. We used a model for linkage analysis that incorporated a high phenocopy rate and a conservative estimate of penetrance. Our goal in this study was not to establish definitive linkage but rather to detect all regions possibly harboring major genes for BP-I in these pedigrees. To facilitate this aim, we evaluated the degree to which markers that were informative in our data set provided coverage of each genome region; we estimate that at least 94% of the genome has been covered, at a predesignated threshold determined through prior linkage simulation analyses. We report here the results of our genome screen for BP-I loci and indicate several regions that merit further study, including segments in 18q, 18p, and 11p, in which suggestive lod scores were observed for two or more contiguous markers. Isolated lod scores that exceeded our thresholds in one or both families also occurred on chromosomes 1, 2, 3, 4, 5, 7, 13, 15, 16, and 17. Interesting regions highlighted in this genome screen will be followed up using linkage disequilibrium (LD) methods.

Plasmodium falciparum antigenic diversity: Evidence of clonal population structure

Plasmodium falciparum, the agent of malignant malaria, is one of mankind’s most severe scourges. Efforts to develop preventive vaccines or remedial drugs are handicapped by the parasite’s rapid evolution of drug resistance and protective antigens. We examine 25 DNA sequences of the gene coding for the highly polymorphic antigenic circumsporozoite protein. We observe total absence of silent nucleotide variation in the two nonrepeated regions of the gene. We propose that this absence reflects a recent origin (within several thousand years) of the world populations of P. falciparum from a single individual; the amino acid polymorphisms observed in these nonrepeat regions would result from strong natural selection. Analysis of these polymorphisms indicates that: (i) the incidence of recombination events does not increase with nucleotide distance; (ii) the strength of linkage disequilibrium between nucleotides is also independent of distance; and (iii) haplotypes in the two nonrepeat regions are correlated with one another, but not with the central repeat region they span. We propose two hypotheses: (i) variation in the highly polymorphic central repeat region arises by mitotic intragenic recombination, and (ii) the population structure of P. falciparum is clonal—a state of affairs that persists in spite of the necessary stage of physiological sexuality that the parasite must sustain in the mosquito vector to complete its life cycle.

Genetic epidemiology of single-nucleotide polymorphisms

On the causal hypothesis, most genetic determinants of disease are single-nucleotide polymorphisms (SNPs) that are likely to be selected as markers for positional cloning. On the proximity hypothesis, most disease determinants will not be included among markers but may be detected through linkage disequilibrium with other SNPs. In that event, allelic association among SNPs is an essential factor in positional cloning. Recent simulation based on monotonic population expansion suggests that useful association does not usually extend beyond 3 kb. This is contradicted by significant disequilibrium at much greater distances, with corresponding reduction in the number of SNPs required for a cost-effective genome scan. A plausible explanation is that cyclical expansions follow population bottlenecks that establish new disequilibria. Data on more than 1,000 locus pairs indicate that most disequilibria trace to the Neolithic, with no apparent difference between haplotypes that are random or selected through a major disease gene. Short duration may be characteristic of alleles contributing to disease susceptibility and haplotypes characteristic of particular ethnic groups. Alleles that are highly polymorphic in all ethnic groups may be older, neutral, or advantageous, in weak disequilibrium with nearby markers, and therefore less useful for positional cloning of disease genes. Significant disequilibrium at large distance makes the number of suitably chosen SNPs required for genome screening as small as 30,000, or 1 per 100 kb, with greater density (including less common SNPs) reserved for candidate regions.

Minisatellite marker analysis of Trypanosoma brucei: Reconciliation of clonal, panmictic, and epidemic population genetic structures

The African trypanosome, Trypanosoma brucei, has been shown to undergo genetic exchange in the laboratory, but controversy exists as to the role of genetic exchange in natural populations. Much of the analysis to date has been derived from isoenzyme or randomly amplified polymorphic DNA data with parasite material from a range of hosts and geographical locations. These markers fail to distinguish between the human infective (T. b. rhodesiense) and nonhuman infective (T. b. brucei) “subspecies” so that parasites derived from hosts other than humans potentially contain both subspecies. To overcome some of the inherent problems with the use of such markers and diverse populations, we have analyzed a well-defined population from a discrete geographical location (Busoga, Uganda) using three recently described minisatellite markers. The parasites were primarily isolated from humans and cattle with the latter isolates further characterized by their ability to resist lysis by human serum (equivalent to human infectivity). The minisatellite markers show high levels of polymorphism, and from the data obtained we conclude that T. b. rhodesiense is genetically isolated from T. b. brucei and can be unambiguously identified by its multilocus genotype. Analysis of the genotype frequencies in the separated T. b. brucei and T. b. rhodesiense populations shows the former has an epidemic population structure whereas the latter is clonal. This finding suggests that the strong linkage disequilibrium observed in previous analyses, where human and nonhuman infective trypanosomes were not distinguished, results from the treatment of two genetically isolated populations as a single population.

DNA variation at the Sod locus of Drosophila melanogaster: An unfolding story of natural selection

Patterns of variation at the Sod locus of Drosophila melanogaster suggest that the protein polymorphism at this locus has very recently arisen. In addition, it appears that a previously rare DNA variant has been recently and rapidly driven to intermediate frequency. From the size of the region (>20 kb) that has been swept along with this rare variant, and patterns of linkage disequilibrium in the region, it is inferred that strength of selection was large (s > 0.01) and that the sweep occurred more than 25,000 generations ago. In addition, there are striking similarities to patterns of variation observed at the Est6 and Est-P loci, which are located approximately 1,000 kb from Sod.

Evidence for clonal propagation in natural isolates of Plasmodium falciparum from Venezuela

We have analyzed 75 isolates of Plasmodium falciparum, collected in Venezuela during both the dry (November) and rainy (May–July) seasons, with a range of genetic markers including antigen genes and 14 random amplified polymorphic DNA (RAPD) primers. Thirteen P. falciparum stocks from Kenya and four other Plasmodium species are included in the analysis for comparison. Cross-hybridization shows that the 14 RAPD primers reveal 14 separate regions of the parasite's genome. The P. falciparum isolates are a monophyletic clade, significantly different from the other Plasmodium species. We identify three RAPD characters that could be useful as “tags” for rapid species identification. The Venezuelan genotypes fall into two discrete genetic subdivisions associated with either the dry or the rainy season; the isolates collected in the rainy season exhibit greater genetic diversity. There is significant linkage disequilibrium in each seasonal subsample and in the full sample. In contrast, no linkage disequilibrium is detected in the African sample. These results support the hypothesis that the population structure of P. falciparum in Venezuela, but not in Africa, is predominantly clonal. However, the impact of genetic recombination on Venezuelan P. falciparum seems higher than in parasitic species with long-term clonal evolution like Trypanosoma cruzi, the agent of Chagas' disease. The genetic structure of the Venezuelan samples is similar to that of Escherichia coli, a bacterium that propagates clonally, with occasional genetic recombination.

Molecular markers reveal that population structure of the human pathogen Candida albicans exhibits both clonality and recombination.

The life history of Candida albicans presents an enigma: this species is thought to be exclusively asexual, yet strains show extensive phenotypic variation. To address the population genetics of C. albicans, we developed a genetic typing method for codominant single-locus markers by screening randomly amplified DNA for single-strand conformation polymorphisms. DNA fragments amplified by arbitrary primers were initially screened for single-strand conformation polymorphisms and later sequenced using locus-specific primers. A total of 12 single base mutations and insertions were detected from six out of eight PCR fragments. Patterns of sequence-level polymorphism observed for individual strains detected considerable heterozygosity at the DNA sequence level, supporting the view that most C. albicans strains are diploid. Population genetic analyses of 52 natural isolates from Duke University Medical Center provide evidence for both clonality and recombination in C. albicans. Evidence for clonality is supported by the presence of several overrepresented genotypes, as well as by deviation of genotypic frequencies from random (Hardy-Weinberg) expectations. However, tests for nonrandom association of alleles across loci reveal less evidence for linkage disequilibrium than expected for strictly clonal populations. Although C. albicans populations are primarily clonal, evidence for recombination suggests that sexual reproduction or some other form of genetic exchange occurs in this species.

An increased risk of Crohn's disease in individuals who inherit the HLA class II DRB3*0301 allele.

The role of inflammatory T cells in Crohn's disease suggests that inherited variations in major histocompatibility complex (MHC) class II genes may be of pathogenetic importance in inflammatory bowel disease. The absence of consistent and strong associations with MHC class II genes in Caucasian patients with inflammatory bowel disease probably reflects the use of less precise typing approaches and the failure to type certain loci by any means. A PCR-sequence-specific oligonucleotide-based approach was used to type individual alleles of the HLA class II DRB1, DRB3, DRB4, and DRB5 loci in 40 patients with ulcerative colitis, 42 Crohn's disease patients, and 93 ethnically matched healthy controls. Detailed molecular typing of the above alleles has previously not been reported in patients with inflammatory bowel disease. A highly significant positive association with the HLA-DRB3*0301 allele was observed in patients with Crohn's disease (P = 0.0004) but not in patients with ulcerative colitis. The relative risk for this association was 7.04. Other less significant HLA class II associations were also noted in patients with Crohn's disease. One of these associations involved the HLA-DRB1*1302 allele, which is known to be in linkage disequilibrium with HLA-DRB3*0301. These data suggest that a single allele of an infrequently typed HLA class II locus is strongly associated with Crohn's disease and that MHC class II molecules may be important in its pathogenesis.

A pre-Columbian Y chromosome-specific transition and its implications for human evolutionary history.

A polymorphic C-->T transition located on the human Y chromosome was found by the systematic comparative sequencing of Y-specific sequence-tagged sites by denaturing high-performance liquid chromatography. The results of genotyping representative global indigenous populations indicate that the locus is polymorphic exclusively within the Western Hemisphere. The pre-Columbian T allele occurs at > 90% frequency within the native South and Central American populations examined, while its occurrence in North America is approximately 50%. Concomitant genotyping at the polymorphic tetranucleotide microsatellite DYS19 locus revealed that the C-->T mutation displayed significant linkage disequilibrium with the 186-bp allele. The data suggest a single origin of linguistically diverse native Americans with subsequent haplotype differentiation within radiating indigenous populations as well as post-Columbian European and African gene flow. The mutation may have originated either in North America at a very early time during the expansion or before it, in the ancestral population(s) from which all Americans may have originated. The analysis of linkage of the DYS199 and the DYS19 tetranucleotide loci suggests that the C-->T mutation may have occurred around 30,000 years ago. We estimate the nucleotide diversity over 4.2 kb of the nonrecombining portion of the Y chromosome to be 0.00014. compared to autosomes, the majority of variation is due to the smaller effective population size of the Y chromosome rather than selective sweeps. There begins to emerge a pattern of pronounced geographical localization of Y-specific nucleotide substitution polymorphisms.

Maintenance of pre-mRNA secondary structure by epistatic selection.

Linkage disequilibrium between polymorphisms in a natural population may result from various evolutionary forces, including random genetic drift due to sampling of gametes during reproduction, restricted migration between subpopulations in a subdivided population, or epistatic selection. In this report, we present evidence that the majority of significant linkage disequilibria observed in introns of the alcohol dehydrogenase locus (Adh) of Drosophila pseudoobscura are due to epistatic selection maintaining secondary structure of precursor mRNA (pre-mRNA). Based on phylogenetic-comparative analysis and a likelihood approach, we propose secondary structure models of Adh pre-mRNA for the regions of the adult intron and intron 2 where clustering of linkage disequilibria has been observed. Furthermore, we applied the likelihood ratio test to the phylogenetically predicted secondary structure in intron 1. In contrast to the other two structures, polymorphisms associated with the more conserved stem-loop structure of intron 1 are in low frequency, and linkage disequilibria have not been observed. These findings are qualitatively consistent with a model of compensatory fitness interactions. This model assumes that mutations disrupting pairing in a secondary structural element are individually deleterious if they destabilize a functionally important structure; a second "compensatory" mutation, however, may restabilize the structure and restore fitness.

Genetic analysis of type 1 diabetes using whole genome approaches.

Whole genome linkage analysis of type 1 diabetes using affected sib pair families and semi-automated genotyping and data capture procedures has shown how type 1 diabetes is inherited. A major proportion of clustering of the disease in families can be accounted for by sharing of alleles at susceptibility loci in the major histocompatibility complex on chromosome 6 (IDDM1) and at a minimum of 11 other loci on nine chromosomes. Primary etiological components of IDDM1, the HLA-DQB1 and -DRB1 class II immune response genes, and of IDDM2, the minisatellite repeat sequence in the 5' regulatory region of the insulin gene on chromosome 11p15, have been identified. Identification of the other loci will involve linkage disequilibrium mapping and sequencing of candidate genes in regions of linkage.

Genome scan of human systemic lupus erythematosus: Evidence for linkage on chromosome 1q in African-American pedigrees

Systemic lupus erythematosus (SLE) is an autoimmune disorder characterized by production of autoantibodies against intracellular antigens including DNA, ribosomal P, Ro (SS-A), La (SS-B), and the spliceosome. Etiology is suspected to involve genetic and environmental factors. Evidence of genetic involvement includes: associations with HLA-DR3, HLA-DR2, Fcγ receptors (FcγR) IIA and IIIA, and hereditary complement component deficiencies, as well as familial aggregation, monozygotic twin concordance >20%, λs > 10, purported linkage at 1q41–42, and inbred mouse strains that consistently develop lupus. We have completed a genome scan in 94 extended multiplex pedigrees by using model-based linkage analysis. Potential [log10 of the odds for linkage (lod) > 2.0] SLE loci have been identified at chromosomes 1q41, 1q23, and 11q14–23 in African-Americans; 14q11, 4p15, 11q25, 2q32, 19q13, 6q26–27, and 12p12–11 in European-Americans; and 1q23, 13q32, 20q13, and 1q31 in all pedigrees combined. An effect for the FcγRIIA candidate polymorphism) at 1q23 (lod = 3.37 in African-Americans) is syntenic with linkage in a murine model of lupus. Sib-pair and multipoint nonparametric analyses also support linkage (P < 0.05) at nine loci detected by using two-point lod score analysis (lod > 2.0). Our results are consistent with the presumed complexity of genetic susceptibility to SLE and illustrate racial origin is likely to influence the specific nature of these genetic effects.